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ABSTRACT 

Aims. We study the dependence of turbulent transport coefficients, such as the components of the a tensor (otij) and the turbulent 
magnetic diffusivity tensor {rjij), on shear and magnetic Reynolds number in the presence of helical forcing. 

Methods. We use three-dimensional direct numerical simulations with periodic boundary conditions and measure the turbulent trans- 
port coefficients using the kinematic test field method. In all cases the magnetic Prandtl number is taken as unity. 
Results. We find that with increasing shear the diagonal components of otij quench, whereas those of rjij increase. The antisymmetric 
parts of both tensors increase with increasing shear. We also propose a simple expression for the turbulent pumping velocity (or j 
effect). This pumping velocity is proportional to the kinetic helicity of the turbulence and the vorticity of the mean flow. For negative 
helicity, i.e. for a positive trace of ajj, it points in the direction of the mean vorticity, i.e. perpendicular to the plane of the shear 
flow. Our simulations support this expression for low shear and magnetic Reynolds number. The transport coefficients depend on the 
wavenumber of the mean flow in a Lorentzian fashion, just as for non-shearing turbulence. 
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1. Introduction 

Understanding the origin of non-uniform large-scale magnetic 
fields in stars, galaxies, and accretion discs continues to pose im- 
portant challenges. Such fields are commonly thought to be the 
result of dynamo action that converts the kinetic energy of tur- 
bulent motions and large-scale shear into magnetic energy. The 
usual framework for s tudying such dynamo act ions is mean field 
electrodynamics (e.g. iKrause & Radlerl Il980b . Over the years, 
however, the applic a bility of this framework has bee n questioned 
dPiddingtonl 1 198 it IVainshtein & Cattaneol Il992l) . In particu- 
lar, an important debate in this c onnection revolves around th e 
role played by magnetic helicity dGruzinov & Diamond!. [1994). 
Magnetic helicity is conserved for ideal (non-dissipative) mag- 
netohydrodynamics (MHD) and also in the limit of large mag- 
netic Reynolds numbers. This conservation places severe con- 
straints on the growth of the mean magnetic field and may reg- 
ulate the quenching of th e a effect as the magnetic Reynold s 
number (Rcm) increases dBrandenburg & S ubramanian. 2005). 
Under certain circumstances (no magnetic helicity fluxes), a 
quenching becomes more severe with a proportional to Re^ 1 . 
It has been suggested that this can be alle viated by expelling 
magn e tic helicity thro u gh ope n boundaries (Blackman & Field, 
2000; Kleeor inet all l2000l). possibly mediated by shear 
dVishniac & Choi 1200 U ISubramanian & Brandenburg! 120041 
iBrandenburg & Sandinl 120041) . Furthermore, shear itself is an 
important ingredient in MHD dynamos in solar and stellar set- 
tings as, for example, in the solar tachocline. Hence, it is im- 
portant to understand how shear affects the turbulent transport 
coefficients, in particular the components of the onj and rjij ten- 
sors. 

Several studies have recently looked at various as- 
pects of thi s problem (see , e.g., [Rogachevskii & Kleeorin, 
120031 120041 [Brandenburg; 120051; iRtidiger & Kitchatrnovl 



2006 ; iRadler & Stepanovi l2006t iLeprovost & Kuril 120081 
Brandenb urg et all l2008al) . These works employ different tools 
and make different assumptions and are often applicable to 
limited regions of the parameter space. As a consequence, care 
must be taken in comparing these results. For example, using 
semi-analytical tools, w hich treat the nonl i near L orentz force 
feedback perturbatively, ILeprovost & Kiml d2008l) have found 
that shear can reduce a in helically forced turbulence. This is 
analogous to the "rotational quenching" of turbulent transport 
coefficients with increasing Coriolis or inver se Rossby numbers 
dKitchatinov et all 1 19941 IPipin et all Il996l) . There are also 
other related cases in which the presence of shear enhances the 
growth rate of the dynamo. For example; using direct numerical 
simulations of the MHD equations in the pre sence of shear 
and non-helical forcing, Youse fetalJ d2008allbb have found 
large-scale dynamos whose growth rate increases linearly with 
shear. Such scaling has also been found for a-shear dynamos 
where the a effect is due to stratified convection with shear 
(Kapy la et all [2008a). Furthermor e, using the kinematic test 
field method (described below), Brandenbur g et all (|2008a) 
have studied the dynamo coefficients in the presence of shear, 
but in the absence of helicity, and they find that Gaussian 
fluctuations of a about zero are strong enough to driye an 
incoherent q- shear dynamo dVishniac & B randenburg. 1997; 
iProctorl [2007). The s ignificance of the incoher ent a effect has 
also been stressed by iHughes & Proctor! ([2008), although their 
system may have also had a net a effect. They dismissed this on 
the grounds that for an imposed uniform magnetic field a is very 
small. However, this result di sagrees with recent calculations of 
a using the test field method (Ka pyla et alll2008bl) . 

In this paper we use three-dimensional direct numerical sim- 
ulations of the kinematic test field equations with helical forcing 
in order to study the dependence of turbulent transport coeffi- 
cients on shear and magnetic Reynolds number. In Section 2 we 
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give a brief account of our model. Section 3 contains our results 
and finally we conclude in Section 4. 



2. The model 



We use the test field method dSchrinner et all I2005L 120071) to 
calculate the turbulent transport coefficients. This method and its 
modification in the presenc e of large-scale shear are described in 
dBrandenburg et all l2008al) . Here we just point out the essence 
of the method and elaborate only on those aspects where our 
treatment differs from their paper. 

In the presence of large-scale shear, the equations of magne- 
tohydrodynamics are treated in the following way. Writing the 
velocity as the sum U + U , where the large-scale shear ve- 
locity U s = Sxy with a constant shear S, and assuming an 
isothermal equation of state characterised by the sound speed c s , 
the momentum equation becomes 



VU 
~Vt 

Here F v 



U VU- SU x y - c 2 s V lnp + / + F vis 



(1) 



V ■ (2pvS) is the viscous force, Sy = 
j V • U is the traceless rate of strain ten- 
sor (not to be confused with the shear parameter S), v is the 
kinematic viscosity, p is the fluid density, c s is the isothermal 
sound speed and 

V d d 

— = — + Sx—. (2) 



Vt dt 



dy' 



As our external forcing / we employ helical, w hite-in-time, ran- 
dom forcing described in ( Brandenburg! l200"lt) . In this paper we 
consider the purely kinematic problem, so there is no feedback 
in Eq. (Q~|l due to the Lorentz force. In addition to the momentum 
equation, we have the continuity equation 



Vp 
Vt 



= -U Vlnp- V-U, 



and the uncurled induction equation in the Weyl gauge, 
VA 

—— = -SAyX + U x B - rjpoJ. 



(3) 



(4) 



Here the magnetic field is B = V x A, the current density 
is J = V x B /no, and no and rj are the vacuum permeabil- 
ity and the molecular magnetic diffusivity, respectively. In the 
mean field approach to MHD one usually decomposes the mag- 
netic field (or magnetic vector potential) and velocity into mean 
(indicated by an overbar) and fluctuating parts respectively 



A = A + a, U = U + u. 



(5) 



The equations satisfied by the mean and fluctuating parts of the 
magnetic vector potential are given by 



VA 
~Vt 



= -SA y x + U x B + £ - i]p J, (6) 



and 

Va 

Vt 



-Sa y x +Uxb+uxB+uxb-£- r]p, j, (7) 



where j = J — J and £ = u x b is the mean turbulent elec- 
tromotive force. Specification of this second order quantity in 



terms of the mean field constitutes a closure problem. A com- 
mon procedure is to expand £ in terms of the mean field B and 
its derivatives, 



£i — £(ii 



(8) 



where a,j and rjijk are the tensorial turbulent transport coeffi- 
cients, and £q quantifies additional contributions that arise even 
in the absence of a mean field, owing to small-scale dynamo ac- 
tion, for example. Here and throughout, summation is assumed 
over repeated indices. In the test field method we take the mean 
magnetic field to be a given 'test field' B, and calculate these 
tensors by measuring £ by solving Eq. (0 simultaneously with 
Eqs. (Q]i and (O, while using Eq. (f5]) to find u. In order to find all 
the components of the ay and r/ij tensors, and not just projec- 
tions relevant to the actual fields, one has to use an orthogonal 
set of test fields and solve Eq. (0 for each of them. 

In our simulations we employ averages over x and y direc- 
tions to define our mean fields, which are therefore functions of 
z and t only. Thus, all other components of Bj^ except B XtZ 
and B ViZ vanish, and the a and rj tensors can be written as rank 
two tensors with rjn = 77^ 3 and r) i2 = -rjm for i,j = (1,2) 
dBran denburg et all l2008al) . Now consider, as an example, the 
following two test fields, 

B~f = B (cos kz, 0,0), (i jf = fcS (0, - sinfcz, 0), (9) 
B i = B (smkz, 0,0), p Ji = kB Q (0, + coskz,0), (10) 

which we use to compute the two corresponding mean electro- 
motive forces £ cl and £ . The relevant components of the a 
and 77 tensors are then given by 



an 
-f]%ik 



cos kz sin kz 
— sin kz cos kz 




(ID 



The other 2+2 components of the a and 77 tensors can be simi- 
larly determined by using the test fields, 

T$f = l?o(0,cosfe,0), B~f = fcfl o (0,sinfo,0). (12) 

In what follows we denote a test field by B Pq , where p = c,s 
and q = 1, 2. A particular small-scale magnetic vector potential 
that develops in response to the B Pq is denoted by a pq and the 
corresponding small-scale magnetic field is given by b pq = V x 
a pq . Note that at large values of ReM there will also be small- 
scale dynamo action that will lead to spurious time-dependencies 
of £. However, since neither the test fields nor or rjij depend 
on time, i.e. 

£f(;z,-t) = £oi{z,t) + a lj B j pq {z) +r) i jkBj i k" I (z) , (13) 

such time dependence must be entirely due to £q 9 (z, t) and can 
be eliminated by time averaging. 

In the following we shall, to begin with, use k = k\, the 
wavenumber corresponding to the box size, to study the depen- 
dence of onj and 77^ on shear and magnetic Reynolds number. 
We shall discuss the dependence of the <x%j and rjij tensors on 
k in Section [331 Note that the usual a pproach of using uniform 
applied fields for calculating a (e.g., Courvois ier et all [2006) 
corresponds to a special case of the test field method for k = 0. 
However, dynamos generate large-scale fields with non-zero fc, 
so it is important to relax this restriction. It is then also important 
to calculate rjij. The test field method allows the simultaneous 
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calculation of all the components of the oiy and rjij tensors for 
arbitrary values of k. 

The test field metho d has recently been criticised by 
Cattaneo & Hughes] d2008l) on the grounds that the test fields are 
arbitrary predetermined mean fields. They argue that the result- 
ing turbulent transport coefficients will only be approximations 
to the true values unless the test fields are close to the actual 
mean fields - a criticism equally applicable to other methods 
using arbitrary uniform fields. Thi s concern has been addressed 
by iTilgner & Brandenburg! d2008l) . who argue that Eq. (0 can 
instead be applied to any mean field. This statement has been 
numerically verified in three cases that we describe below. 

Firstly, the test field method correctly r eproduces a vanish- 
ing gr owth rate in saturated nonlinear cases dBrandenburg et all 
2008b). Secondly, in the time-dependent case, the test field 
method correctly reproduces also a non-vanishing growth rate, 
but in that case it is no longer permissible to express £ in 
terms of a multiplication of turbulent transport coefficients with 
the mean field and its spatial derivat ives. One must therefore 
write Eq. ((H) as a convolution in time dHubbard & Brand"en burg, 
2008). Finally, the success of the test field method becomes 
particularly clear when it is applied to a passive vector field 
that obeys a separate induct ion equation with a velocity field 
from a saturated dynamo (Tilgne r & Brandenburg! 120081). T his 
question was originally posed by ICattaneo & Tobiasl d2008l) . In 
particular for the Roberts flow with a mean field of Beltrami 
type, e.g. one that is proportional to (cos kiz, sinfciz, 0), the 
a,ij tensor is anisotropic and has an additional component pro- 
portional to BiBj that tends to quench the components of the 
isotropic part of CKy. The fastest growing passive vector field 
is then proportional to (sin kxz, — cos kiz, 0). This result has 
been c onfirmed both numerically and using weakly nonlinear 
theory dTilgner & Brandenburg! 120081) . 

In the following we ignore the complications involving time- 
dependent mean fields and restrict ourselves to transport coeffi- 
cients that apply strictly speaking only to the time-independent 
or marginally excited case. For our numerical simulations we use 
the Pencil Cod^J], where the test field algorithm has already 
been implemented. All our numerical simulations are performed 
in a periodic cubic box. The forcing scale is chosen to have the 
wavenumber kf/ki = 5. This gives enough scale separation for 
a large-scale field to develop (Haugen et al. 2004), and is still 
not too big to reduce the resulting Reynolds numbers too much. 
We use units such that c s = k\ = po = Mo = 1 and arrange the 
forcing amplitude so that the Mach number is around 0.1. For 
the magnetic Prandtl number we choose Ptm = vjr\ = 1, where 
v lies in the range 3 x 10~ 3 to 2 x 10 -4 (in units of c s /k\). We 
choose the shear S such that the parameter 

Sh = S/{u rms k { ) (14) 

takes values in the range —0.02 to —0.9, where u lms is the root- 
mean-square velocity. All the runs are started with uniform den- 
sity, p = po, U — 0, and a pq — 0. Depending upon the pa- 
rameters of a particular run we use up to 256 3 grid points. For 
each test field B we calculate the time averages of the a>ij and 
rjij tensors over time intervals over which it rms is statistically 
stationary. In runs with higher magnetic Reynolds numbers, 

Re M = Wrms/(??&f), (15) 

we obtain an exponential growth of small-scale magnetic field 
(see below). We interpret this as being associated with the 

1 |http : / / www . nordita .org/ software /pencil- code| 



£ p ) q (z,t) term. This often gives rise to large fluctuations in all 
components of the a« and rjij tensors at late times. In such 
cases we confine our calculations of the time averages to time 
intervals over which 6^ s does not exceed B Pq by more than a 
factor of about 20. Up until this point, the components of the 
ctij and rjij tensors show only their intrinsic fluctuations, but at 
later times these will be swamped by additional contributions 
from the small-scale dynamo that grows exponentially in time. 
If necessary, we repeat our calculations over several independent 
realizations by resetting a pq = at regular time intervals. 

3. Results 

Our principal aim in this paper is to study the effects of varying 
shear (Sh) and magnetic Reynolds number (Rcm) on the com- 
ponents of the Qy and rjij tensors. In the subsection below we 
summarize our results concerning the different components of 
these tensors. 

3. 1. Diagonal components of the transport tensors 

The isotropic parts of the o>ij and rjij tensors are respectively 
characterised as 

a = Kan + a 22 ), lit = 5(17x1 + V22}, (16) 

where (•) denotes an average taken over z and t. We normalize 
these quantities by = — |it rms and ?/ t o = ^u lms /kf which 
are their respective expressions obtained using the First Order 
Smoothing Approximation (FOSA) for ReM <C 1, which has 
previously been confirmed with the test field method in s imula- 
tions of helical turbulence without shear dSur et all 12008). 

Since u Tms enters the normalization of <Xy and rjij, it is use- 
ful to first look at how it changes as a function of ReM (Fig-QJ) 
and Sh (Fig. QJ). As can be seen, u rms increases as a function 
of ReM for small ReM ~ 1 and then reaches a plateau for high 
ReM- On the other hand, u rms is almost a constant as a function 
of shear except for high Sh, where we observe excitation of the 
vorticity dynamo discussed further in Section [3~4l 

In Fig.|2]we show a as a function of shear for three different 
values of the magnetic Reynolds number ReM ~ 1, 20 and 72. 
Figure [2}} shows the corresponding results for rj t . We note that 
for small shear the turbulent transport coefficients are close to 
their values for zero shear. As |Sh| increases, a decreases and rjt 
increases up to four times rjto- As can be seen there is a clear 
reduction (quenching) of a with increasing shear in all these 
cases. In order to examine the possible convergence of the re- 
sults with ReM we plot a and rjt as functions of ReM for three 
different values of the shear parameter Sh, see Fig. [3] For small 
ReM ~ 1 we observe an increase in both a and rj t with ReM- A 
similar initial increase of a and rjt was also seen in earlier simu- 
lations of helical turbulence without shear dSur et all 120081) an d 
in non-helical shear flow turbulence (Brandenburg et al., 2008a). 
For higher values of ReM and Sh, both a and rjt show large vari- 
ations. In these kinematic simulations, however, we e xpect them 
to ten d to constant values asymptotically at high ReM dSur et all 
120081) . 

In an e arlier study of shear flow t urbulence with non-helical 
forcing by Brandenburg et al. (2008a), the diagonal components 
of rjij were found to be the same. Obviously, in the absence of 
helicity all the components of are zero. It turns out that in 
the presence of helicity the two diagonal components of rjij are 
still the same, but those of onj are now non-zero and not equal 
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Fig. 1. Root-mean-square velocity u rms normalized by the speed of sound c s as: (a) a function of ReM for different values of Sh: 
Sh ~ 0.07(*), 0.2(A), andO. 3(D) respectively and (b) a function of the shear parameter Sh for different value of ReM: ReM ~ 1 
(*), Re M ~ 20 (A), and Re M ~ 72 (□). 




Fig. 2. Turbulent transport coefficients (a) a/a and (b) r/t/^to as functions of shear parameter Sh for different values of ReM: 
Re M ~ 1 (*)> ReM ~ 20 (A), and Re M ~ 72 (□). 




Fig.3. (a) ajoiQ and (b) ?7t/??to as functions of ReM, for different values of Sh: Sh ~ 0.07(*), 0.2(A), andO. 3(d) respectively. 
Horizontal dashed lines at a/a = 1 and fft/ffto = 1 are added to facilitate comparison. 
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to each other. This is best shown by considering the quantities 

( a = |(an - "22), e,, = 5(1711 ~ 1722}- (17) 

The results are shown in Fig. H] Note that especially for large 
values of ReM the values of e a are predominantly negative. Since 
both an and 022 are negative, this means that \an | is larger than 
|o;22 1, although the relative difference is only of the order of at 
most 5 per cent. 

3.2. Off-diagonal components 

We now consider the off-diagonal components of atij and rjij . 
The results are depicted in Figs. [5] and [6] Of particular in- 
terest among these is the component 7721, whose dependence 
on Sh and ReM is shown in Figs. [5Jl and [6}l, respectively. 
This component can indicate the possible presence of a shear- 
current dynamo that may operate whe n r)2iS/(r]Tki) 2 > 1 
dRogachevskii & Kleeorinl I2003L |2004|) . Here r] T = T] t + r) is 
the total (sum of turbulent and microscopic) magnetic diffusiv- 
ity. In our case we have S < 0, which implies that the necessary 
condition for the shear-current dynamo to operate is 7721 < 0. 
As can be seen from Fig. [6ll, for the range of ReM values con- 
sidered here 7721 is positive for small shear but becomes nega- 
tive for strong shear and ce rtain values of ReM- Earlier work of 
[Brandenburg et al. (2008a) without helicity did indicate a sim- 
ilar sign change, although only for larger ReM- However, the 
error bars were so large that this result could not be regarded 
as significant. For the run with the strongest shear (— Sh « 0.3) 
and for ReM = 40 we now find 7721 to be more clearly negative, 
but for smaller and larger values of ReM the results are again, at 
least within error bars, compatible with zero. Also, of course, the 
present results apply to the case wi th helicity and are therefor e 
not really comparable with those of Brandenburg et al. (2008a), 
where the helicity is zero. The antisymmetric contributions to 
the ccy and 77^ tensors are characterised by the vectors 

7fc = -\tijkUij, 5 k = -^eijkTjij. (18) 

Since our averages depend only on z, the z components of these 
tensors are irrelevant and therefore only the z components of the 
7 and 8 vectors are of interest. We denote those simply by 7 and 
6, with 

7 = 2"( a 2i ~ CC12), £=5(1721-7712). (19) 

A time series of both quantities is shown in Fig. [7] for positive 
and negative signs of the kinetic helicity. The results show that, 
in our case with S < 0, 7 is positive (negative) for positive 
(negative) kinetic helicity, whilst 8 has always the same sign. 

Some idea about the functional forms of 7 and 8 can be ob- 
tained from symmetry considerations. The vectors 7 and 6 enter 
the electromotive force thus 

£ = ...+7XB — fi 6 x J, (20) 

so we see that 7 must be a polar vector and 8 must be an axial 

g 

vector. Using the shear flow U , the only axial vector that can be 

5 

constructed is the mean vorticity, W = V x U , so we expect 8 
to have a component that is proportional to W. Likewise, since 
7 is a polar vector which points in the direction of W, the two 
must be related via a pseudoscalar. In the present case the only 
pseudoscalar available is the kinetic helicity, u> ■ u. Based on 
these symmetry arguments we write 

7 = CyT 2 LJ-uW, 8 = C S T 2 rfW : (21) 



where we have introduced two non-dimensional quantities, C 7 
and Cs which could be either positive or negative, and r is a 
correlation time that we approximate here by 

T = (Urmsfcf)" 1 . (22) 

In the present case, u> ■ u itself is a negative multiple of a. We 
therefore expect 7 to have a component proportional to aW, 
multiplied by the correlation time r. By similar arguments we 
expect S to have a component proportional to %W, multiplied 
by t. Based on these arguments we can also write 

7 = C 7 oWt, 8 = C s %Wt, (23) 

with new non-dimensional quantities, C 1 w C 7 /3 and C$ « 
C5/3, that we expect to be of order unity. Our simulations con- 
firm this reasoning and show that both coefficients are of order 
unity with positive C 1 and negative Cs (e.g., C 1 — 0.5 and 
Cs = 0.25 for the run shown in Fig. [8]) for runs with small ReM 
and Sh; see Figs. [8] and [9l respectively. The sign in Eq. [231 has 
been chosen a posteriori from our simulations. For high shear 
and Reynolds number this simple reasoning is no longer accu- 
rate and, at least in one case (Fig |9^), C 7 even becomes nega- 
tive. For lower values of Sh we find an almost constant negative 
7 of the order of 0.05ao (or less). This result is fairly indepen- 
dent of ReM- The moduli of 7 and S increase with increasing 
shear parameter, and for larger values of ReM both quantities 
may approach an asymptotic value. For small ReM we observe 
a roughly linear increase of 7 with shear (see inset of Fig [8^), 
further verifying Eq. ( |2T1 >. 

The coefficient 7 can be interpreted as turbulent pumping, 
i.e. advection of the magnetic field by means other than the 
mean velocity field. In strongly stratified convection, turbulent 
pumping has been seen to be d i rected from higher to lower tur- 
bulence intensi ty (Tobia s~et al.L [1998, 2001; Oss endriiver et all 
2002; Ka pyla et al.L 120061) . which is usually in the downward 
direction. Thus, turbulent pumping is likely to play an impor- 
tant role in convection zones of stars where it can overcome the 
buoyancy of the magnetic field. In the present case where the tur- 
bulence is homogeneous, however, stratification does not play a 
role and the pumping is just due to the combined action of shear 
and helical turbulence. 

3.3. Scale-dependence 

So far we have confined the calculations of the ctij and 77^ ten- 
sors to test fields whose characteristic length scale is the largest 
scale in the domain, i.e., test fields of the form sin kz or cos kz 
with k = ki, where fci = 2ir/L and L is the box size of our sim- 
ulations. We have also done similar calculations for other values 
of k, which correspond to test fields with different characteristic 
length scales. In Fig.[l0]we show the dependence of a and rj t on 
k for ReM ~ 3. The decrease can be modelled by a Lorentzian 
peaked at k = (Fig. [TOt. This result is not surprising, because 
Loren tzian fits have been obtaine d earlier also in the absence of 
shear dBrandenburg et all l2008d) . 

3.4. Effects from the vorticity dynamo 

In the runs with relatively high values of shear and Reynolds 
number, the root-mean-square velocity, U lms , which initially 
reaches a steady state, shows an exponential growth at late times, 
as can be seen in Fig. [TTa . Similar behaviour is seen for tur- 
bulent transport coefficients, which also show large fluctuations 
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Fig. 4. Normalized (a) e a and (b) e v ), as defined in Eq. d!7t . as functions of ReM for constant value of the shear parameter Sh: 
Sh ~ 0.07(*), 0.2(A), and 0.3(D) respectively. Horizontal dashed lines are added to facilitate comparison. 




0.01 0.10 1.00 0.01 0.10 1.00 




Fig. 5. Normalized off-diagonal components of the a and i] tensors as functions of the shear parameter Sh for different values of 
ReM: ReM ~ 1 (*), ReM ~ 20 (A), and ReM ~ 72 (□). Horizontal dashed lines are added to facilitate comparison. 



for later times. Large fluctuations of the tu rbulent t r ansport co- 
efficients at later times are also seen by ISur et all (120081) . but 
those fluctuations are more irregular and have a different ori- 
gin. They are interpreted as being due to the development of a 
small-scale dynamo at late times and are not due to the more 
systematic increase in the rms velocity, which in turn is associ- 
ated with the shear. However, we emphasize that even before b pq 
becomes dominated by this type of dynamo action, the temporal 
and spatial fluctuations of ay and rjij are of the order of oq and 
77tn, respectively. This is tr ue even for large values of ReM; see 
Brandenburg et al. (2008a), who have argued that these fluctua- 
tions also contribute to large-scale dynamo action via the inco- 



herent a effect dVishniac & Brandenburg! 1 1997tlProctor[ 12007). 
For example, the time-series of r\% shown in the inset of Fig.[TT1 
has a plateau even beyond the range over which C/ rms is steady. 

The vertical spikes in the inset of Fig. QT|come from reset- 
ting a' pq to zero in regular time intervals; see the end of Sect. [2] 
Note, however, that the mean values and also the upper enve- 
lope trace the evolution of U lms reasonably well, including the 
increased rise after t — 300. The late time behaviour is accompa- 
nied by the formation of large-scale vortical structures, as seen 
in Fig. [TTp. which is a signa ture of th e vorticity dynamo p ro- 
posed bv lElperin etail d2003l) . see also dYousef et all l2008allbl) . 
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Fig. 6. Normalized off-diagonal components of the a and r\ tensors as functions of ReM for constant values of Sh: Sh ~ 
0.07 (*), 0.2 (A), and 0.3 (□). Horizontal dashed lines are added to facilitate comparison. 



A detailed numerical study of the vorticity dynamo has been per- 
formed by Kap vla et al. (200l. 

The presence of the vorticity dynamo and the resulting sys- 
tematic variation in U lms as well as the turbulent transport co- 
efficients for late times limit the lengths of time over which av- 
erage values of the components of the ctij and r)ij tensors can 
be calculated. In the above case this interval lies between the 
two arrows in Fig [TTk This limits the range of Sh and ReM 
that we have probed and explains the larger error bars shown 
for example in Fig. [6] and why they cannot be reduced by 
simply running our simulations for longer times. This problem 
would be avoided in the presence of magnetic fields, because 
then th e resulting Lorentz force woul d quench the vorticity dy- 
namo (Kapvla & Brandenburg, 2008). This is however beyond 
the scope of this paper. 

4. Conclusions 

We have studied the effects of varying shear and magnetic 
Reynolds number on the turbulent transport coefficients in the 
presence of helicity in the kinematic limit. We have shown that 
for fixed ReM, a is reduced (quenched) with increasing shear. 
Despite the differences in the assumptions made, this quench- 
ing is qualitative l y sim ilar to the recent results obtained by 
iLeprovost & Kiml |2008). To the best of our knowledge this is 
the first numerical study to show quenching of a as a func- 
tion of shear in helical turbulence. We find that r/ t increases 
with increasing shear in the range of ReM values that we have 
considered here. A simil ar behaviour for rjt was also seen in 
[Brandenburg et al. (2008a) where the forcing was non-helical. 

We also compute the off-diagonal components of the a,j 
and ifij tensors. The antisymmetric part of a.^ corresponds to a 



turbulent pumping velocity 7 in the direction perpendicular to 
the plane of the shear flow. It shows a roughly linear increase 
with shear for small magnetic Reynolds number. We propose 
simple expressions for 7 and S in Eq. (f2Tb which show reason- 
able agreement with our numerical results for small Sh and ReM- 
Our expression shows that for negative helicity, 7 points in the 
direction of the vorticity of the mean flow. Regarding the com- 
ponent 7721 we find indications that, at least in one or two cases, 
this component changes sign and becomes negative. This could 
be of significance in connection with the shear-current effect. 

We also find that all the turbulent transport coeffi- 
cients depend on the wavenumber of the mean flow in a 
Lorent zian fashion, just as in the case of non-shearing turbu- 
lence (Brandenburg et al., 2008c). This means that the kine- 
matic values of a and rj t for k = k\ are close to the 
values obtained for k — -> 0. This is not the case for cer- 
tain non-turbulent flows such as the Gallowa y-Proctor flow 
(ICourvoisierl [2 008: Radl er & Brandenburg! 120081) . In an earlier 
paper. lCourvoisier et al.ld2006l) considered only the limiting case 
k = for this flow. 

Several aspects of the present investigations could be of as- 
trophysical relevance. Turbulence in celestial bodies is helical 
and exhibits an a effect. In addition, shear (S) can be an im- 
portant ingredient in that the efficiency of large-scale dynamo 
action is determined by the product of a and S. However, 
as S increases, a itself becomes quenched when S becomes 
comparable with the inverse turnover time, i.e. St = 0(1). 
Furthermore, the turbulent diffusivity becomes enhanced, sup- 
pressing the dynamo even further. Finally, it is found that the 
combined action of helicity and shear gives rise to a pumping 
velocity of mean magnetic field perpendicular to the plane of the 
shear flow. The existence of such a pumping velocity has not 
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Fig. 7. Normalized time series of turbulent transport coefficients, (a) a/|ao|, (b) rj t /r]to, (c) 7/|«o| and (d) S/r]to, from two runs 
with exactly the same parameters (ReM ~ 1.4, Sh = —0.5), but different signs of helicity for the external force. Continuous and 
broken lines denote results from runs with positive and negative helicities, respectively. 




Fig. 8. Normalized (a) 7 and (b) 6 as functions of Sh for ReM ~ 1 (*), ReM ~ 20 (A), and ReM ~ 72 (□). Horizontal dashed lines 
are added to facilitate comparison. The inset shows 7/|ao| versus — Sh for ReM ~ 1. 



been emphasized before. On the other hand, it is well known 
that a £l (or rather aS) dynamo s can have travelling wave solu- 
tions dBrandenburg et all |2001|) . When the product of a and S 
is positive, these waves travel in the positive z direction, which 
agrees with the direction of pumping. In the near-surface shear 



layer of the Sun this pumping would therefore support the equa- 
torward migration in that layer. 

It is important to understand the quenching of turbulent 
transport coefficients in the presence of shear beyond the kine- 
matic approximation. In that case one needs to include the 
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Fig. 9. Normalized (a) 7 and (b) 5 as functions of Rejyr for constant values of the shear parameter Sh: Sh ~ 
0.07 (*), 0.2 (A), and 0.3 (□). Horizontal dashed lines are added to facilitate comparison. 




Fig. 10. (a) a(k)/ao and (b) rjt(k)/r]to as a function of the wavenumber k of the large-scale magnetic field for ReM ~ 3. The solid 
lines show the Lorentzian fits. 



induction equation (in addition to the test field equations) 
and incorpora te the resulting Lorentz force in the momen- 
tum equation (Brandenburg et al., 2008b). This would also help 
in alleviating problems of strong late-time fluctuations aris- 
ing from the vorticity dynamo, because the vorticity dynamo 
tends to be suppressed by magnetic fields of equipartition 
strengths (Kapyla & Brandenburg, 120081) . Similarly, given their 
potential importance in allowing the escape of magnetic helicity, 
the effects of open boundary conditions also needs to be con- 
sidered. These questions are under study and will be reported 
elsewhere. 
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